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Abstract. The dynamical system approach has recently acquired great 
importance in the investigation on higher order theories of gravity. In this talk I 
review the main results and I give brief comments on the perspectives for further 
developments. 
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1. Introduction 

Since the very first proposal of General Relativity alternative formulations of the law of 
gravity have been proposed for various purposes: from unification of the fundamental 
interaction to the explanation of the dark energy phenomenon. In spite of the great 
variety of proposals, all these attempts present a common drawback: they are much 
more complicated than General Relativity. This means that it is very difficult to devise 
tests on these theories, and hence there are not many opportunities to gain insights 
on the their physics. 

Obtaining such information is critical to the debate about the nature of 
gravitation and its relation with the dark energy issue. Particularly important in 
these respect are the indirect methods, i.e. methods that allows to solve indirectly 
the gravitational field equations or to obtain a qualitative idea of how these solutions 
behave. 

One very interesting methods of this type is based on the application of the 
dynamical system approach to cosmology. This method has been developed for general 
relativity since the 70s, and it has lead to beautiful insights into the evolution of the 
anisotropic cosmologies pQ. 

In this paper we summarize the results obtained from the application of dynamical 
system theory to fourth order gravity and in particular to i?™-gravity. We also discuss 
very briefly some further applications of this method. 

Unless otherwise specified, we will use natural units (H — c — ks = 8ttG = 1) 
throughout the paper and Greek indices run from to 3. The semicolon represents 
the usual covariant derivative and the "dot" corresponds to time differentiation. 



2. R n — gravity 



The action for the gravitational interaction in this theory reads 

A = J d^x^-g [ X (n)R n + L M ] , (1) 

where x( n ) is a positive function of n that reduces to 1 for n — 1. For R ^ 0, the field 
equations for this theory can be written as 
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where T^f, is the stress energy tensor for the standard matter. In this way, the non- 
Einstcinian part of the gravitational interaction can be modelled as an effective fluid 
which, in general, presents thermodynamic properties different from standard matter 
U In the following we will analyze this equation in the Friedmann-Lemaitre-Robertson- 
Walker (FLRW) and Bianchi I metrics with the aid of the dynamical system approach. 

X Note that the nature of the matter term and the indctcrmination in the sign of R makes this theory 
to be fully meaningful only if we consider in the following only values of n belonging to the set of the 
relative numbers Z and the subset of the rational numbers Q, which can be expressed as fractions 
with an odd denominator. 
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3. Dynamical analysis of the FLRW case 

In the FLRW metric, the ([2]) take the form 



with 
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(i + 3ifyt(l + w) = . (5) 



i? = - 6 (- + i + 4) ( 6 ) 



where = a/ a, n is the spatial curvature index and we have assumed standard matter 
to be a perfect fluid with a barotropic index w. Note that in the above equations we 
have considered a and R as two independent fields so that the equation are of effective 
order two and the conservation equation for matter is the same as standard GR [2]. 

The form of the above equations suggests the following choice of expansion 
normalized variables: 

x = m (n - lh y = 6^ (1 ~ n) ' z = ^wW-i ■ x = (7) 

Differentiating with respect to the logarithmic time Af — In a, we obtain the 
autonomous system (matter is considered a perfect fluid) 



2(n — 2)y 9 xy 
- 2a; - 2a; 2 y — + (1 + a;)z - 3zw 



n — 1 n — 1 

y' = — ^— [(3 - 2n)x - 2y + 2(n - l)z + 2(n - 1)1 . (8) 
re — 1 

z' = z ^2z - 1 - 3a; ^ - 3w^ ,l + a;-y-z + ; ; s: = 0. 

Note that the two planes y = and z = correspond to two invariant submanifolds. 
This implies that for this system no finite global attractor exists. The behavior of the 
scale factor corresponding the fixed points can be found using the equation 

H=(x c + - z c - lj H 2 . (9) 

If n 1 and [n — l)(a;c — zc — 1) — yc 7^ this equation can be integrated to give 

(Uc \ 
1 — xc h zc) ■ (10) 
n-1 J 

In the same variables the energy density /x can be written as 

/i = zy"- 1 ^ 2 ™ , (11) 

thus for n > 1 both the [x, y) and [x, z) planes are invariant vacuum manifolds, but 
if n < 1 the vacuum submanifold is not necessarily compact. 

A detailed analysis of this dynamical system is shown elsewhere [5] . Here we will 
focus on the interval 1.367 < n < 2. This interval is suggested by a fitting of the data 
coming from WMAP and Supernovae la data [3]. 
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Table 1. Coordinates of the finite and asymptotic fixed points and solutions for 
FLRW /{"-gravity in vacuum. 



Point 



Coordinates 



Solution/Behaviour 



A 
B 
C 
V 



[x,y] = [0,0] 
[x,y] = [-1,0] 

\ T „1 - \ 2{n-2) 4n-5 
J/J [ 2n-l ' 2n— 1 

[x,y] = [2(1 -n) ,2(n-l) 2 ] 
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3.0.1. The vacuum Case In a vacuum spacetime (i.e. /i = 0), the variable z is 
identically zero and the third equation of ([5]) becomes an identity. Let us analyze this 
case first. Setting x' = and y' = 0, we obtain the four fixed points shown in Table 
[U In our interval the solution associate with the point C is a power law inflation and 
the only finite attractor. This is an interesting result because this model contains 
cosmic histories that naturally approach to a phase of accelerate expansion. However, 
the presence of the invariant submanifolds already mentioned makes this attractor not 
global and we have to check if there are other attractors in asymptotic regime. 

The idea that the system above might have a nontrivial asymptotic structure 
comes from the fact that the phase plane is not compact. The asymptotic analysis 
can be easily performed using the Poincare approach [4] . We obtained six fixed points 
which are summarized with their behavior in Table [TJ The stability analysis shows 
that in the interval of n we have chosen there is another attractor (this time global): 
the point T> x . Using the asymptotic limit of Equation[9](see [2]), we found that D M is 
associated with a Lemaitre type evolution in which the universe reaches a maximum 
size and then recollapses. The presence of "D^ reduces the measure of the set of 
initial conditions for which an orbit will approach to C, or in simpler words it makes 
this type of cosmic histories "less probable" . A pictorial representation of the whole 
(compactified) phase space is given in Figure [TJ 

3.0.2. The non-vacuum Case When matter is present we consider the full system 
©. Setting x' = 0, y' = 0, z' = and solving for x,y,z we obtain seven fixed 
points and two fixed subspaces. The first four points lay on the plane z — and 
correspond to the vacuum fixed points. Of the other three fixed points, Q is definitely 
the most interesting because of its associated solution: a = a t 3(1 + m > , fi — ^t . 
The presence of such a fixed point brings the idea that there could be cosmic histories 
in this model in which a(n unstable) Friedmann-likc phase is followed naturally by a 
phase of accelerated expansion. Applying the standard tools of the dynamical system, 
we discover that this is actually the case in our interval for n. In fact, for 1.367 < n < 2 
the point C is an attractor and the point Q is a saddle, they are both physical and 
placed in a connected sector of the phase space. Thus in principle a cosmic history like 
the one pictured above is possible. This is also confirmed by numerical investigation 
of the dynamical system. The last thing to check is the presence of other attractors 
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in the asymptotic regimes. The compactification is achieved in the same way as the 
vacuum case, but the equations obtained are much more complicated. Here we will 
limit ourselves to say that there are other attractors in the phase space, in the form of 
subspaces and single points and their presence reduce the measure of the set of initial 
conditions that lead to a cosmic history connecting Q and C. 



4. Dynamical analysis of the LRS Bianchi I cosmologies 



In order to obtain the simplest possible form for the field equations in the LRS Bianchi 
I metric, we use the 1+3 covariant approach to cosmology [5]. In this formalism the 
cosmological equations can be written as 
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(12) 
(13) 
(14) 



together with the ([5]). Here is the volume expansion = 3^ = ZH and a is the 
square root of the magnitude the symmetric shear tensor er a {,. The set of expansion 
normalized variables is: 

3i? , 3-R . 3m „ 3ct 2 



RO y " 11 a 2n8 2V xnR n - x Q 2 
and the dynamical system can be written as 
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(16) 



(2 - Zw) - z + S 



n - 1 



1 — Y, + x — y — 2 = 0. 
The solution associated to every fixed point can be found via the equations 



e = 



n- 1 



yi - - 2 



e 2 



a = -- (2 + Ej + + Zj) 9 a. 



(17) 



For n^l, (7^0 and a R.H.S. different from zero, these equations admit the solution 

-l 



a = a (t- t )" 



a = (Too ' 







'2 + E, -^t//,) . (18) 

-(2 + E, +y 4 + z,) . (19) 

As in the previous section we will focus on a specific intervals of n referring the reader 
to [6] for a general analysis. 
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Table 2. The finite an asymptotic fixed points and eigenvalues for ij n -gravity in 
a LRS Bianchi I model. 
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4-0.3. The vacuum Case In order to treat the vacuum case, we set z = and we 
neglect the equation for z as in the FLRW case. Setting £' = and y' = we 
obtain the fixed subspaces and their solutions (see Table HJ. Note that the fixed point 
A represents the same solution associated to the point A. Direct verification with 
the cosmological equations reveals that this point is anisotropic i.e. a a b = 0. The 
presence of such a point is particularly interesting because it means that there could 
be cosmic histories in which there is an isotropic state for an otherwise anisotropic 
universe. Two interesting cases arise. When 1 < n < 5/4, the point A is a repeller 
and part of the orbits emerging from it approach the fixed line. This means that 
in these cosmic histories the universe starts in an isotropic state and then develops 
anisotropies (see figure [1} . The presence of an isotropic past attractor is not present 
in GR and implies that in i? n -gravity, like in brane cosmology, there is no need for 
special initial conditions for inflation to begin [6]. When n > 5/4, A is an (isotropic) 
attractor and all the orbits above the fixed line start on C\ and converge to it (see 
Figure [T]) . This scenario presents the same smooth transition between a decelerated 
and accelerated expansion of the FLRW case. If we analyze the asymptotic regime for 
this last case we can see that for n > 5/4 the point is an attractor (see Table [5]), 
but is in a separate section of the phase space with respect to A, sot that A can be 
considered an "effective" global attractor. 

4-0.4- The Matter Case Setting £' = 0, y' = and z' = in (fH)]) we obtain, together 
with the vacuum fixed points, other two fixed points, one of which (Q) is an isotropic 
point associated with the same solution of Q. Let us check if it is possible to have 
a cosmic history similar to the one we found in FLRW. For n > 5/4, the fixed line 
contains repulsive points nearby the origin, the point Q is unstable and the point A 
in an attractor. Since these points are not separated by any invariant subspace, there 
is in principle an orbit that connects them. Along this orbit, the universe would start 
in an anisotropic state, evolve towards a more isotropic state to smoothly approach to 
a phase of accelerated expansion. Of course the issue of determining how many other 
attractors are present (and so how "probable" this evolution is) is still present and 
we have to perform a detailed asymptotic analysis in order to understand this point. 
Using the results of [B] it is possible to show that there are other attractors in the 
phase space that might influence the global evolution. However, the outcome is that 
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Figure 1. Global (compactified) phase spaces for B n -gravity. The first picture 
represent the FLRW case with 1.36 < n < 2; the second the Bianchi I case for 
1 < n < 5/4 and the third the Bianchi In > 5/4 case. The standard capital 
letters represent the finite fixed points the underlined capital letters represent the 
asymptotic fixed points. The point E is nor a fixed point but it is defined as the 
last attractive point of the fixed line. The curly capital "1" represents the line Cl- 



an orbit like the one described above is still possible and surely deserves more study. 
5. Conclusion 

In this paper we hope to have shown how useful and powerful the dynamical system 
approach is in dealing with complicated model like fourth order gravity ones. The 
application of this method has allowed new insights on higher order cosmological 
models and has shown a deep connection between these theories and the cosmic 
acceleration phenomenon, which is worth to be further studied. 

The perspective for future application of the dynamical system approach to higher 
order gravity can be basically divided in two thrusts. The first one is to analyze more 
complicated Lagrangians. Some of work in this direction has already started in [7] with 
quadratic gravity. The second one is the generalization to more complicated metrics, 
which will allow us to consider different physical framework. In the future, both of 
these thrusts will make possible not only to develop experimental test for alternative 
gravity but also to allow a better understanding of the reasons underlying the success 
of General Relativity. 
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